Optimized Monotonic Convex Pair Potentials Stabilize Low-Coordinated Crystals 
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We have previously used inverse statistical-mechanical methods to optimize isotropic pair interac- 
tions with multiple extrema to yield low-coordinated crystal classical ground states (e.g., honeycomb 
and diamond structures) in d-dimensional Euclidean space R d . Here we demonstrate the counterin- 
tuitive result that no extrema are required to produce such low-coordinated classical ground states. 
Specifically, we show that monotonic convex pair potentials can be optimized to yield classical 
ground states that are the square and honeycomb crystals in R 2 over a non-zero number density 
range. Such interactions may be feasible to achieve experimentally using colloids and polymers. 



The forward approach of statistical mechanics focuses 
on finding the structure and macroscopic properties of 
many-particle systems with specified interactions. This 
approach has led to the discovery of rich and complex 
many-particle configurations The power of the in- 
verse statistical-mechanical approach is that it can be 
employed to design interactions that yield a targeted 
many-particle configuration with desirable bulk physical 
properties.— This work continues our general program to 
use inverse approaches to optimize pair interactions to 
achieve novel targeted classical ground-state configura- 
tions in d-dimensional Euclidean space R d . In particular, 
we have found optimized pair interactions that yield low- 
coordinated crystal classical ground states (e.g., square 
and honeycomb crystals 5 - in R 2 , and simple cubio^ and 
diamond^ crystals in R 3 ), materials with negative ther- 
mal expansion, 8 negative Poisson's ratio^ and designed 
optical properties. lp We envision using colloids and/or 
polymers to realize such designed potentials because one 
can tune their interactions &Hr— 

Earlier uses of the inverse approach 5 - did not regard 
experimental feasibility as a constraint. These investiga- 
tions allowed a largely unconstrained class of spherically 
symmetric pair potentials. In some instances in which 
the goal was to target low-coordinated crystal ground 
states, it was shown that only a few potential wells were 
required^ which nonetheless may be difficult to realize 
experimentally. If purely repulsive monotonic pair poten- 
tials existed that could achieve unusual ground states, 
they would be easier to produce experimentally. How- 
ever, encoding information in monotonic potentials to 
yield low-coordinated ground-state configurations in Eu- 
clidean spaces is highly nontrivial. Such potentials must 
not only avoid close-packed (highly coordinated) com- 
petitors but crystal configurations that are infinitesimally 
close in structure (very slight deformations of the tar- 
geted low-coordinated crystal), which is a great challenge 
to achieve theoretically while maintaining the monotonic- 
ity property. 

In this Letter, we use a modified inverse approach 
to obtain monotonic convex potentials whose ground 
states in R 2 are either the square lattice or honeycomb 
crystal?^ Thus, our work is a theoretical proof of con- 



cept that monotonic convex potentials can stabilize low- 
coordinated crystals. 1 We begin by describing briefly 
the procedure that we employ to optimize the mono- 
tonic convex potentials for the targeted low-coordinated 
crystals aided by generalized coordination functions; see 
Rcf. 15. This is followed by an analysis of their stability 
characteristics. 

We consider the total potential energy <&n(t n ) of a 
configuration C of N particles with positions r 1? r 2 , ... 
to be given by a sum of pairwise terms: 

M^HXX^)' (1) 

where v(r) is the isotropic pair potential and = |r$ — 
Yj | . For a targeted configuration C* to be a ground state 
associated with a potential v, the total potential energy 
needs to satisfy the following property: 

*n(v,C*)<$ n (v,C) VC. (2) 

By expressing the potential v as a function of M parame- 
ters ai, aM, i.e., v = v(ai, . . . , am), it should be possi- 
ble to find the optimized potential by varying the param- 
eters until inequality @ is satisfied for all possible con- 
figurations C. Due the uncountably infinite number of 
possible configurations, it is impossible to check them all. 
Instead, we restrict ourselves to a subset of them, which 
we call the competitor configurations C. This allows us 
to redefine the problem as an optimization, 5 where the 
objective is to maximize the energy difference between 
the targeted configuration C* and its closest competi- 
tors. This is done by introducing a utility variable A 



We consider monotonic potentials that are also convex because 
this class of interactions may be amenable to rigorous analysis. 
Indeed, in Ref. 14, Cohn and Kumar have rigorously constructed 
potentials that stabilize unusual targeted configurations on the 
surface of a d-dimensional sphere using only monotonic convex 
pair potentials. Restriction to compact spaces made their prob- 
lem much easier to solve because their pair potentials had com- 
pact support set by the sphere radius. Nonetheless, their results 
are suggestive that similar proofs can be constructed in R d . 
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which is to be maximized while satisfying the following 
constraints: 

^ N (v( ai ,...,a M ),C*) < $at (v( ai ,...,a M ),C)-A VC € 

(3) 

For a fixed potential v, the utility variable A can only 
be as large as the smallest energy difference, $n {v, C) — 
<&at (v, C*), between a competitor and the targeted con- 
figuration. Since the functional form v is allowed to vary, 
the optimization procedure will find the potential that 
maximizes the energy difference between the targeted 
configuration and its closest competitor. 

For a given targeted configuration C* , we begin with a 
competitors set C that only includes the triangular lat- 
tice. A trial pair potential v\ is optimized using that set, 
and its putative ground state Ci is computed using the 
Metropolis Monte Carlo algorithm. Namely, we attempt 
to determine the ground state for the trial potential v\ by 
generating an initial configuration from a Poisson point 
process. This configuration is then slowly annealed using 
the Metropolis scheme down to zero temperature. Since 
this algorithm cannot guarantee that the obtained con- 
figuration is the ground state, we repeat the procedure 
multiple times and keep the lowest-energy configuration 
as the trial ground state C\ of the trial potential v\ . 

If the energy of G\ is lower than that of C* for the 
trial potential, it proves that C* is not the ground state 
of vi . To discriminate against potentials with C\ as their 
ground states, we add C\ to the list of competitors C, be- 
fore optimizing a new trial potential vi. This procedure 
is repeated until we arc confident that the trial-potential 
ground state is indeed C* , at which point we have found 
our optimized potential. This method is adapted from 
the one presented by Cohn and Kumar with the dif- 
ference being that we only add to C configurations that 
have lower energies than C* for a given trial potential. 

We restrict ourselves to v(r) that are sums of M = 12 
negative powers of r, with a cutoff at r = R > 0: 



M 



v(r) 



i=l 

r> R. 



(4) 



Additionally, we only consider continuous potentials 
whose first and second derivatives are also continuous at 
the cutoff. These two conditions guarantee that the inter- 
action forces — V r Nii are continuous and that the phonon 
spectra can be calculated, respectively. Furthermore, a 
scale is imposed by setting v[r = 1) = 1. 

The first low-coordinated crystal configuration to 
be targeted with our simulated-annealing optimization 
method is the square lattice with a nearest-neighbor 
distance of unity subject to the condition that <2j s 
[—1000,1000]. The number density for such a configu- 
ration is unity (p = 1). The optimization procedure is 
restricted to monotonic convex pair potentials that are 
zero beyond a cutoff distance r = R = 2. We find the 




FIG. 1: (Color online) (a 



Optimized convex pair potential 
targeting the square lattice. The potential v(r) at unit dis- 
tance (where r is measured in terms of the nearest-neighbor 
distance) is taken to be unity. The vertical lines represent the 
nearest-neighbor distances for the honeycomb crystal (black 
solid), the square lattice (red dashed) and triangular lattice 
(blue dotted) at a number density of unity (p = 1). |(b)| The 
optimized potential ground state, obtained by slowly anneal- 
ing the systems, starting from a fluid. The annealing was 
performed in a 20 x 20 box containing 400 particles under 
periodic boundary conditions. For illustration purposes, the 
point particles are shown to have finite sizes. 



following optimized pair potential: 



v(r) = < 
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r < 2, 

r > 2, 
(5) 

To confirm that the ground 
state of potential ([5]) is indeed the square lattice, we per- 
formed multiple simulated-annealing calculations. All of 
them resulted in either square lattices (as shown in Fig. 
1(b)) or slightly deformed square lattices, whose ener- 



which is plotted in Fig. 1(a) 



gies were always higher than that of the perfect square 
lattice. We use M = 12 terms in our potential because 
higher M cause numerical instabilities and lower M re- 
sult in potentials that only weakly discriminate against 
competitors. Importantly, the potential function ([S]) is 
only one example within a large class of functions that 
could be optimized to stabilize the square lattice. 

To see intuitively why the purely repulsive potential 
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(b) 



FIG. 2: (Color online) As in Fig. [T] except that it is for the 



optimized potential targeting the honeycomb crystal. In (a) 
the ne arest -neighbor distances are calculated for p — 4/3 \73 
and in |(b)| we use 416 particles in a periodic box with dimen- 
sions 24 x 13^3. 



(O succeeds in stabilizing the square lattice consider the 
interactions due to the first and/or second coordination 
shells, the main contributors to the total energy. For the 
square, triangular and honeycomb crystals (with p = 1), 
the first coordination shell contributions to twice the 
total energies per particle u [u = Yljjti v ( r ij) f° r tne 
studied crystals] are respectively given by 4 x v(l) = 4, 
6 x u(1.075) = 4.762 and 3 x w(0.877) = 4.963 (these 
can be compared to the entire function it, which is re- 
spectively 4.456, 4.764 and 5.236). We see that the lower 
value of the optimized potential ([5]) for the triangular lat- 
tice at its nearest-neighbor distance is not enough to com- 
pensate for the higher coordination number compared 
to the square lattice (six instead of only four). Thus, 
the relatively slow decrease of the potential §5§ around 
r = 1 allows it to favor configurations with low coordi- 
nation, even if it means having closer nearest neighbors. 
The lower-coordinated honeycomb crystal is also discrim- 
inated against, due the large increase of v(r) for r < 1. 

However, there are more subtle configurations that 
have to be discriminated against beside the aforemen- 
tioned ones. These include infinitesimally close con- 
figurations, such as the rhombical and rectangular lat- 
tices. The difference in u between a rectangular lat- 
tice of aspect ratio 1 + e and the square lattice is 
equal to [«'(l)/2 + s/2v'(V2) + v"{l)/2]e 2 + 0(e 3 ) = 
2.551e 2 -f 0(e 3 ). The difference in u between a rhom- 
bical lattice of angle ir/2 — e and the square lattice is 



[v'(l)+v'(V2)/V2+v"(V2)}e 2 + 0(e 3 ) = 2.110e 2 +O(e 3 ). 
Stabilizing the square lattice against these two very close 
neighboring configurations is thus an equilibrium be- 
tween having large second derivatives at the two first co- 
ordination shells, while having preventing the first deriva- 
tive to be too negative at these two shells. 

The second targeted ground-state configuration is the 
honeycomb crystal with a nearest-neighbor distance of 
unity and number density p = 4/3v3- This is a more 
challenging ground state to achieve with a monotonic 
convex potential because it is only trivalently coordi- 
nated. The optimization procedure is still restricted to 
monotonic convex pair potentials, but the cutoff is set to 
R = 3. The optimized pair potential is given by 
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r < 3, 

r > 3, 
(6) 

which is plotted in Fig. |2(a) As for potential (JSJ) , we 
confirmed that the ground state of potential ([6]) is indeed 
the honeycomb crystal by performing multiple simulated- 
annealing calculations. Figure 2(b) shows the result of 



one of those runs that converged to the honeycomb crys- 
tal. All final configurations other than the honeycomb 
crystal had higher energies than that of the perfect hon- 
eycomb, which is strong numerical evidence that the hon- 
eycomb is indeed the ground state of potential ©• 

As in the square-lattice case, the ability of the purely 
repulsive potential ((6]) to stabilize the low-coordinated 
honeycomb crystal lies in its slow decrease near r = 1. 
Consequently, the contributions to u from the first coor- 
dination shells of the honeycomb, triangular and square 
crystals (at p — 4/3\/3), which are respectively 3, 4.077 
and 3.196. The low coordination of the honeycomb crys- 
tal thus compensate for closer neighbors. The second 
shell energy contributions of the square and triangular 
lattices turn out to be larger than that for the honey- 
comb crystal, even if all them are relatively small due 
to the rapid decrease of v(r). We have also verified that 
potential ([5]) discriminates against slightly sheared de- 
formations of the honeycomb crystal, which is consistent 
with our phonon analysis below. 

We also studied the phonon characteristics of poten- 
tials ([5]) and ([S]), i.e., the mechanical response of the 



crystals to small deformations. Figures 3(a) and 3(b) 
show the squared frequency of phonon modes as a func- 
tion of their wave vectors for the square and honeycomb 
crystals. The absence of any negative squared frequency 
indicates that all of the modes have real frequency; thus, 
the crystals are stable under small deformations. 

We have further studied the stability of our optimized 
potentials by exploring the effects of adding point defects 
to the crystals. Since defects cost energy, our targeted 
ground states are stable under such local modifications. 
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FIG. 3: (Color online) Phonon spectra for the optimized po- 
tentials. The squared phonon frequency uj 2 is plotted in term 
of a representative subset of phonon wave vectors q. Note 
that the spectra were calculated over the entire Brillouin zone, 



and no modes with imaginary frequencies were found, (a) 



Phonon spectrum for the "square-lattice" potential shown in 



Fig. 1(a) (b) Phonon spectrum for the "honeycomb-crystal" 
potential shown in Fig. 2(a) 



We also used the newly introduced generalized coordi- 
nation functions^ to show that potentials (J5J) and © 
are part of a large class of monotonic convex pair po- 
tentials that stabilize the square and honeycomb crystals 
and thus our potentials are robust against shape change. 
These details are given in our companion paper 

Are the low-coordinated crystal ground states stable 
over a density range around the density values for which 
they are designed (p = 1 for the square lattice and 
p = 4/3\/3 for the honeycomb crystal)? We have per- 
formed both simulated- annealing ground-state calcula- 
tions and computed the phonon spectra at various den- 



sities. We find that for potential (|5|), the square lattice 
is the ground state for the density range p £ [0.96, 1.10]. 
For potential ((6|), the honeycomb crystal is the ground 
state for p £ [0.74,0.80]. Therefore, both targeted con- 
figurations can be stabilized over a non-zero number den- 
sity range, which is a desirable feature for experimental 
realizations of our optimized potentials. This property 
is not at all obvious for system under positive pressure, 
as their are no a priori reasons why deformed lattices 
that maintain a constant nearest neighbor distance, such 
as the rhomboidal and rectangular lattices, are not the 
ground states for densities other than the one at which 
the optimization was conducted. 

To summarize, whether potentials exist that stabilize 
low-coordinated crystal ground states in Euclidean space 
without any potential wells is not at all obvious. We 
have shown that potentials without wells, namely, mono- 
tonic convex repulsive pair interactions, can produce low- 
coordinated ground states in K 2 , such as the square lat- 
tice and honeycomb crystals. Therefore, the naive ex- 
pectation that purely repulsive interactions will only lead 
to densely packed configurations in Euclidean space has 
been disproved. 

Lindcnblatt et al<^ have fabricated so-called "hairy 
colloids" . These colloids are formed by grafting polymer 
chains onto the surface of nanoscopic microgel spheres, 
in a matrix of polymer chains. Manipulation of such sys- 
tems offer the possibility of mimicking the interactions 
that stabilize the square and honeycomb crystals defined 
by Eqs. §5§ and ([6]), respectively, although such experi- 
mental realizability remains an open fascinating question. 
Note that experimental feasibility should not require con- 
vexity. For example, we have also shown numerically 
that the square lattice can be stabilized by a monotonic 
non convex potential^ In future research, we intend to 
determine whether convex potentials can stabilize three- 
dimensional low-coordinated structures, such as simple 
cubic or diamond crystals. 
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